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Abstract 

We consider the discrete surface growth process with relaxation to the minimum [F. Family, 
J. Phys. A 19 L441, (1986).] as a possible synchronization mechanism on scale-free networks, 
characterized by a degree distribution P{k) ~ k^^, where k is the degree of a node and A his 
broadness, and compare it with the usually applied Edward- Wilkinson process(EW) [S. F. Edwards 
and D. R. Wilkinson, Proc. R. Soc. London Ser. A 381,17 (1982) ]. In spite of both processes 
belong to the same universality class for Euclidean lattices, in this work we demonstrate that for 
scale-free networks with exponents A < 3 the scaling behavior of the roughness in the saturation 
cannot be explained by the EW process. Moreover, we show that for these ubiquitous cases 
the Edward- Wilkinson process enhances spontaneously the synchronization when the system size 
is increased. This non-physical result is mainly due to finite size effects due to the underlying 
network. Contrarily, the discrete surface growth process do not present this flaw and is applicable 
for every A. 
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The study of the dynamics on complex networks became a subject of great interest in 
the last few years since it was realized that they are useful tools to understand biological, 
social and communications systems . Networks are constituted by nodes associated to 
individuals, organizations or computers and by links representing their interactions. The 
classical model for random networks is the Erdos-Renyi (ER) model characterized 
by a Poisson degree distribution P{k) = exp [—(/;;)] {k)^ /k\ where k is the degree or number 
of links that a node has and (k) is the average degree. However, it was found [l] that many 
real networks are characterized by a scale-free (SF) degree distribution given by 

P{k) = (A - 1) /, , (1) 

for kmin < k < kmax whcrc kmin IS the smaller degree that a node can have, kmax is a cutoff 
that diverge when the system size N oo and A represent the broadness of the distribution. 
Most of observed networks such as Internet, the World Wide Web and metabolic networks 
have 2 < A < 3 

It was shown that the topology of the network is very relevant to determine their statics 
and dynamics properties, such as robustness and percolation thresholds , the average 

shortest path length 8] and transport (9|. An important quantity characterizing networks 
is its diameter (maximal hopping) d. In a network of a total of nodes, d scales as IniV 

n □ 

[5|], which leads to the concept of "small worlds" or "six degrees of separation" . For 



scale-free (SF) networks with X < 3 1m d scales as InlnA^, which leads to the concept of 

nn 

ultra small worlds [2|, 112[ . 

Very recently, the research focus is changing from the study of the network topology to 
the study of dynamical processes on the underlying network. Of particular interest are the 
studies on the dynamics and fluctuations of task completion landscapes of queuing networks. 
If for each node on the network there is a scalar h which specifies the time it takes to finish a 
job or the amount of work that has been assigned to it, the fluctuations on h indicates how 
synchronized or balanced is the system. Jobs synchronization and load balance are required 
in many applications such as packet routing on the Internet [l3] or in parallel computing 

These synchronization processes are usually mapped into a non-equilibrium sur: 
growth via an Edwards- Wilkinson (EW) equation IGj on complex networks 17 



ace 



19|. 



The EW equation for the evolution of the growing interface in complex networks is given by 
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—l = y^ Aij{hj -hi)+r]i, (2) 

3=1 

where hi = hi{t) is the height of the interface of node i, Aij is the element ij of the 
adjacency matrix that take the value 1 if z and j are connected and zero otherwise, is the 
system size, z/ is a coefficient that represent the "surface tension" and rji = rii{t) is a random 
Gaussian uncorrelated noise with {r]i} = 0; {f]if]j} = 2DSijS{t — t'), where D is the diffusion 
coefficient and {.} represent averages over configurations. The interface is characterized by 
its roughness W{t) at time t, 

1/2 

Wit) ' ' " 



J.J2(h^~{h)r\ , (3) 

i=l ) 



N 

that represents the fluctuations of the height of the interface around his mean value {h). 

There are several technical advantages of using the continuous EW equation to model 
queue synchronization or load balance processes [l7|, [l8|, [l9| mainly because it is a linear 
continuous equation. However, some real implementations of this processes are intrinsically 
discrete. For this reason, in this work we use the discrete growth model of surface relaxation 



to the minimum (SRM), which is very well known on Euclidean lattices |20|, |2l|, on SF 
networks. It is also well known the fact that on Euclidean lattices this discrete model belongs 
to the same universality class of the EW equation. This might be one of the motivations of 
using EW to model these discrete process. In the SRM model at each time step a node i is 
chosen with probability If we denote by Vi the nodes nearest neighbors of i, then 



if 



hi < hj Vj G f j =^ hi = hi + 1, else, 

node j has the minimum height & Vi =^ hj = hj + 1. 

(4) 



This rules mimics a process where the higher loaded node distributes the excess of load 
to one of his neighbors which is less charged. To generate SF graphs of size A^, we employ 



the Molloy-Reed algorithm (MR) 221]: initially the degree of each node is chosen according 
to a SF distribution, where each node is given a number of open links or "stubs" according 
to its degree. Then, stubs from all nodes of the network are interconnected randomly to 
each other with the two constraints that there are no multiple edges between two nodes and 
that there are no looped edges with identical ends. 



We use for the simulation kmin = 2 because when kmin > 1 there is a high probabihty that 



the network is fully connected 



12| which is required in this work to analyze the interface. 



At t = we initialize all the values of hi with random numbers taken from an uniform 
distribution in [0, 1]. At each time step we select a node with probability and use the 
rules given by Eq. (jl]), then the time is increased by We compute W{t) for SF networks 
with A > 2 and different values of A^. 

In Fig. □ (a) and (b) we plot W{t) for the SRM as function of t for A = 2.5 and A = 3.5 
respectively. In both figures we can see a very short growing regime for W{t) after which the 
system saturates with a width Wg. This fast regime before the saturation can be explained 
in terms of finite size effects. For almost all growth processes the correlation length grows 
with time until it reaches the characteristic length of the system {2^, which for complex 
random networks is the diameter d. As explained above, the diameter is very small, and the 
system reaches the saturation time very fast. We focus the attention on the steady state of 
Ws because only at the steady state matters to analyze the fluctuations in the load balance 
of multiprocessors in parallel computing or synchronization of queues. 

For A = 2.5 we found by a linear fitting of W{t) in the steady state that Wg behaves 
with A^ as Wg ~ In A^ (see the inset of Fig. [11(a)). The same scaling behavior was obtained 
for all other values of A < 3. We also run all the simulations for an initial fiat interface and 
found no differences in Wg {2^. For A = 3.5, Wg does depend weakly on the system size for 
big enough networks (see the inset of the Fig. [11(b)). Korriiss reported this lack of finite size 
effect for the growing network model of Barabasi- Albert [1| that has A = 3 24l |. 

As mentioned above, it is well known that this model in Euclidean lattices belongs to the 
EW universality class represented by Eq. ([2D, so it is expected that Eq. ([2|) will show the 
same scaling behavior as the SRM model. In this work we demonstrate that surprisingly 
this is not generally true. In Fig. [2l (a) and (b) we show W{t) as function of t for different 
values of A^ from the numerical integration of Eq. ([2]) with u = 1 and D = 1 ioi SF networks 
with kmin = 2 for A = 2.5 and A = 3.5. 

Counterintuitive, for A = 2.5 Wg decreases with the system size, which is a non expected 
result for any growth model. If this were the case, increasing the system size will be a simple 
strategy to minimize the roughness and thus improving synchronization of queues or balance 
in the load of multiprocessors in parallel computing. 

We next show that the decreasing of the width for A < 3 is mainly due to finite size 
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effects introduced by the MR construction. It was shown in 2J] that for the EW process in 
unweighted networks the absolute lower bound of the is 



WLn = (1 - • (5) 

The decreasing on the width observed in our numerical results is because (k) increases 
with A^. As a consequence of the MR construction which introduces the natural cutoff 

k = k ■ Ari/(A-i) 



25|, {k{N)) is given by 



1 _ 1 //v(A-2)/(A-l) 

{k{N)) = k^ L__ , (6) 

where k^o = k{N oo). Taking into account the results presented in Eq ([5]) where we 
replace < A; > by < A;(A^) > we propose that 

~ W^(oo) — + - ^ (7) 

vv^y^) ^ ^ jY(A-2)/(A-i) J ' y'J 

where Ws{oo) = Ws{N~^oo). 

The fitting of with Eq.Q shows an excellent agreement [see the inset of Figl2](a)] 
with the simulations supporting that the decrease in the width for A < 3 is mainly due to 
the MR construction and for large A^, ~ cte. 

Thus, the scaling behavior of Ws for the SRM model with A < 3 is not well represented 
by the EW equation with constant coefficients u and D, despite the fact that it is often used 
in synchronization problems. 

Next we analyze finite size effects for A > 3. For the SRM model Wg was well fitted by Eq. 
([7]) [see the inset in Fig. [11(b)]. Thus, in this regime, the finite size effects can be attributed to 
the MR construction. For the EW equation we find the best fitting ~ (oo)(l — A/N) 
with B ^ in Eq 171 This behavior cannot be explained as finite size effects due the MR 
construction because for A > 3, A^ diverges faster than A^(-^~2)/(a-i) may be is due to 
the EW process. The separation between the finite size effect of the EW process and the 
MR construction is still an open question that goes beyond the aim of this paper and could 
be the subject of future researches. 

In summary, we simulate the SRM model in SF networks and compare the results with 
the EW process. We show that a discrete model and a continuous model which share the 
same scaling properties on Euclidean lattices does not exhibit this equivalence on complex 
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networks. For the SRM model in SF networks Wg diverges with the system size as In N for 
A < 3. For A > 3 for both, the model and the EW equation, when N ^ oo, Wg ^ cte. 
In order to compare the results of the SRM model with a continuous equation further 
investigation including higher order of the Laplacian in the continuous equation are needed. 
Also the dynamics could introduce some weigths on the links on the underlying unweighted 
network that even at a linear approximation could affect the EW unweighted process. This 
is the aim of our future research. Finally, we can conclude that despite the fact that the SRM 
model and the EW equation belongs to the same universahty class in Euclidean networks, 
in SF networks they do not have the same behaviour. 
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FIG. 1: Plots of W{t) for SRM and different system size N,N = 256 (Q), N = 512 (□), N = 1024 
(o), iV = 2048 (A), iV = 4096 (<i) and N = 8192 (y) for : (a) A = 2.5 we can see that Ws increase 
with the system size A^. In the inset figure we show Ws as function of N in hnear log scale (O)- 
The dashed lines represent the logarithmic fitting supporting that Wg ~ In A^. (b) A = 3.5 we can 
see that W{t) depend weakly on N. In the inset figure we show in symbols Wg as function of A^. 
The dashed lines represent the fitting of with Eq.d?]) {A ~ 10 and B ~ 0.25). In all the inset 
of data's figures we do not display the errors bars because they are of the size of the symbols. 
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FIG. 2: Plots of W{t) from the integration of the EW equation and differents system size A^.(a) For 
A = 2.5, N = 128 (O), N = 256 (□), N = 512 (o), N = 1024 (A), N = 2048 (y), N = 4096 (x) 
and N = 8192 (*). We can see that Ws decreases with N. In the inset figure we show in symbols 
as function of N. The dashed hue is the fitting with the Eq.© {A ^ 0.10 and B « 0.75). (b) 
For A = 3.5, N = 128 (Q), N = 256 (□), N = 512 (o), N = 1024 (A). The dashed line represent 
the fitting with Eq.©. {A ^ 1.15 and B ^ 0) 
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